YCTP-P5-92 



hepth@xxx/9202001 



(N 
> 



Fractional quantum Hall effect 

and 

nonabelian statistics 



ON 
On 

: N. Read 

CO 



Departments of Applied Physics and Physics, P.O. Box 2157 
Yale University, New Haven, CT 06520 



O . 

^ ' and 
CN 

O 

CN 

5' G. Moore 

Department of Physics, P. 0. Box 6666 
■ Yale University, New Haven, CT 06511 

> 
X 

It is argued that fractional quantum Hall effect wavefunctions can be interpreted as conformal 
blocks of two-dimensional conformal field theory. Fractional statistics can be extended to 
nonabelian statistics and examples can be constructed from conformal field theory. The 
Pfaffian state is related to the 2D Ising model and possesses fractionally charged excitations 
which are predicted to obey nonabelian statistics. 

Talk presented by N. Read at the 4th Yukawa International Symposium, Kyoto, 
Japan, 28 July to 3 August, 1991. 

1 



I. Introduction: Particle Statistics and Conformal Field Theory 

This paper is a brief overview of some aspects of the relationship between the theories of the 
fractional quantum Hall effect (FQHE) [[Ij] and two-dimensional conformal fields (CFT) [[|, 
which has been explored in more detail in || . The present section describes general aspects 
of particle statistics, especially nonabelian statistics; the second section uses the Laughlin 
wavefunctions as an example of the relation with CFT; the third section presents a case of 
nonabelian statistics in a FQHE system (the Pfaffian state). 

The notion of particle statistics in quantum mechanics usually refers to the action of the 
permutation group S n on the wavefunction for a collection of n identical, indistinguishable 
particles: the wavefunction is taken to transform as a definite representation of this group, 
which interchanges particles. The usual examples are Bose statistics and Fermi statistics, 
which are the trivial and alternating representations, respectively. A more modern approach 
prefers to exchange particles along some definite paths in space, and allows for the possibil- 
ity that the wavefunction be not well-defined when particles coincide, so that intersecting 
exchange paths must be omitted. Topology then enters the subject, and while it turns out 
for spatial dimension d > 2 that the group of exchanges still reduces to S n , for d = 2 the 
group of topologically distinct exchanges is Artin's braid group B n 0]. The infinite group 
B n can be generated by a set of elementary exchanges Bi, i = 1, . . . n — 1 which simply 
exchange particles i, i + 1 along a path not enclosing any other particles. The most familiar 
representation of B n is fractional statistics, where each Bi acts on the wavefunction as e t6 
(8 real) and so is a one- dimensional, abelian representation; this includes Bose and Fermi 
statistics as special cases. 

We will use the term "nonabelian statistics" when particle wavefunctions transform as a 
nonabelian representation of the permutation or (especially) of the braid group, and we have 
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introduced the term "nonabelions" for such particles ||. In this case, not all the representa- 
tives of the Bi commute, and must be matrices acting on vector wavef unctions. Thus, even 
when the positions and quantum numbers of the particles have been specified, the wave- 
function is not unique but is a member of a vector space (a subspace of the Hilbert space). 
The vector space may be finite dimensional but its dimension will grow with the number of 
particles. Physical observables are invariant under exchanges, so cannot distinguish which 
state in this space the system is in, just as the phase of the wavefunction cannot be directly 
measured in the usual (one dimensional) case. In the latter case, differences or changes 
of phase, and interference effects, can, however, be observed, and analogously, nonabelian 
statistics can manifest itself in physical effects. 

For the permutation group, nonabelian statistics is an old idea known as "parastatistics" . 
It appears, however, that parastatistics cannot be realized in a local theory in d > 2 in a 
nontrivial way (the literature containing this result is reviewed in 0). The issue is whether 
the global degeneracy and nonlocal effects of exchanges are compatible with the natural 
physical idea of locality of interactions, as realized for example in quantum field theory 
governed by a local Lagrangian density. 

An analogous analysis for the braid group in d = 2 has been carried out in [[J . It involves 
not only the matrices Bi but also the notion of "fusion rules" . A theory will usually contain 
several distinct types of particles, each with its own B matrices, and taking a particle of one 
type around one of another type will produce a matrix effect like that of B 2 . Fusion of two or 
more particles into a composite produces either annihilation of the particles or else a particle 
of some new type with its own statistics properties. The matrix elements for these processes, 
denoted F, describe the fusion rules and will have to satisfy consistency relations involving 
B J s. For example, taking a third particle around a pair before or after they fuse should give 
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the same result, since the third particle is far from the fusing pair and so cannot distinguish 
the close pair from their composite (locality). Frohlich et al || conclude that nonabelian 
statistics can be acceptable in two dimensions. Indeed, local actions are known that produce 
nonabelian properties: they are Chern-Simons terms for nonabelian gauge fields, and the 
particle worldlines are represented as Wilson lines ||. Note that as in the abelian case, the 
statistics effects can be transferred between the wavefunctions and the Hamiltonian by a 
singular gauge transformation; the general discussion above was for the gauge choice where 
statistics is exhibited in the wavefunction, i.e. no Chern-Simons-type vector potentials in 
the Hamiltonian. 

As a partial example of nonabelian statistics, we present the following system. In the 
example, we may create particles a, but only in even numbers. The a particles have no 
internal quantum numbers. For two particles, the vector space for fixed positions Zi,Z2 is 
one dimensional, and the wavefunction is (z\ — 22) 8 , so one might imagine that we have 
abelian fractional statistics 6/tt = —1/8. However, for four particles, the vector space is 
two-dimensional. If we place particles at 0, 1 and oo (which in this example can always be 
done through a global conformal transformation), the wavefunctions are functions only of 
the remaining coordinate z, and (choosing a basis) are 

^ ± (z) = (z(l - z))- 1 / 8 ^l±VT^. (1.1) 

In the z-plane, there are branch points at z = 0, 1. As z is analytically continued around 
z — 1, we see that ip± are interchanged. This operation corresponds to two exchanges of the 
particle at z with that at 1, i.e. to B 2 . Other double exchanges involving z are diagonal 
in this basis. To exhibit B itself we would have to perform a conformal transformation to 
obtain z i— > 1, 1 i— > z. Thus we have a single multibranched function whose two sheets 
are linearly independent functions of z, and so the a particles behave as nonabelions. For 
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general number 2n of a particles the wavefunctions form a 2" _1 -dimensional space. 

The example arises from a particular conformal field theory (CFT), the two dimensional 
Ising model at its critical point. We next describe some general features of CFT ]/]]. In a 
1 + 1-dimensional system with short range interactions (or a local lagrangian) at a critical 
point, one has not only scale invariance but also conformal invariance, an infinite-parameter 
group under which correlation functions transform covariantly. Exactly at the critical point 
(so that corrections due to slowly decaying irrelevant or marginal operators are omitted) the 
correlation functions of a collection of fields {0j r (x r )} can be split in the form (taking the 
"diagonal" case for simplicity): 

(n^W) =El jF p;n...»„( z ir--4)| 2 (1-2) 

\r=l / V 

where z r = x r + iy r . The conformal block functions T v are multibranched functions, analytic 
in their arguments, except possibly when two z's coincide. The variable p labelling different 
functions runs over a finite set in a rational CFT. As the z's are varied so as to exchange some 
0i's, the functions T v are analytically continued to different sheets, but can be expressed as 
z independent linear combinations of the original functions through some braiding matrices 
B as we saw above. In this way, the correlation function can be single valued. 

Another operation that can be performed on the correlation functions is the operator 
product expansion. As the arguments z\, z 2 of two fields approach one another, the operators 
merge into a linear combination of single operators: 

&(z)<P 3 (w)~Y,C*(z-w)Mw) (1-3) 

k 

as z — > w, where <pk is some new field of type k and C^(z — w) is a singular coefficient 
function. This operation can be used to define some new matrices F, the fusion matrix that 
describes which fields k appear in the product of % and j. In || the consistency conditions 
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that must be satisfied by B, F were analyzed; these same conditions emerged later in the 
work of Frohlich et al. Thus CFT provides many examples of nonabelian statistics, when 
conformal blocks are interpreted as wavefuctions. In the following we show that this idea 
applies directly to FQHE states. 

II. Example: Laughlin States 

In this section we show that both the fractional statistics of quasiparticles in the Laughlin 
states U and the actual wavefunctions are related to CFT conformal blocks. The construc- 
tion given here is from but parts of it have been obtained independently by others. 
In 1 + 1-dimensional Euclidean spacetime, define a free scalar field by its correlator, 

{y{z)y{z')) = -\og{z - z') (2.1) 

where the log is complex, so the field creates right-moving excitations only. All correlators 
can be obtained from this one using Wick's theorem. Now consider the function 




where p = l/2irq, q is an integer, and the integral is taken over a disk of area 2irqN centered 
at the origin. We expand and contract using fl2.1|) . Each exponential is assumed normal 
ordered, i.e. y?'s from the same exponential are not to be contracted together. The result is 


Y[( Zi -^j'e-sE,/^ 1 *-'). (2.3) 

i<j 

In the last factor, the real part of the log gives e~* ^ ' 2l ' 2 for Z{ inside the disk. Apart 
from the remaining phase, we now recognize the function as Laughlin's wavefunction for 
particles in the lowest Landau level at filling factor v = 1/q. Since the vertex operators 
expiyfqcp(z) are associated with Coulomb charges, it is clear that we have a holomorphic 
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version of Laughlin's 2D plasma picture ||, which gives the wavefunction and not just its 
modulus squared! 

The imaginary part of the log in ([2.3|) winds by 27rd 2 z' as each Zi goes round each point 



z' in the integration region, so is highly singular [fTI |. It contributes a pure phase to the 
function, which can be gauged away by an equally singular gauge transformation. This 
phase describes a uniform magnetic field, which we identify as the physical magnetic field in 
the FQHE problem. A nice way to see this is to consider not the phase itself but the change 
as one particle describes a closed loop C. Clearly the exponent changes by 



1 

2^ 



d 2 z'2m (2.4) 



where the integral is over the area enclosed by C, so the phase change is 1/27T times the area 
of the loop, which is just the Berry phase for adiabatic transport of a particle of charge 1 in 
a field of strength 1/2tv. In all the following equations, the singular phase will be implicitly 
gauged away. 

In a similar way, one can obtain quasihole wavefunctions as (for two quasiholes): 
/e^^e^^ JJ e Vw>(*0 e -V9 J d 2 z'p<p(z')\ 

= (z- w) 1 ' 11 ]J(z - z k )(w - z k ) Y[(z t - Zj y e-i^ MMMMM 3 . ( 2 .5) 

k i<j 

The first factor gives explicitly the fractional statistics, 9/rr = 1/q. Apparently, the construc- 
tions from CFT correlators give results in the gauge where all Berry phases — background 
magnetic field as well as fractional statistics — appear in the wavefunction, not as a vector 
potential (connection) in the Hamiltonian. This will be important to us in the next section 
when we move on to nonabelian statistics (adiabatic transport). 

We now return to the 2 + 1-dimensional world of electrons in two dimensions in a strong 



magnetic field and describe briefly the order parameter picture of the FQHE, following [12 



(see also [H| |Xlf). In words, a composite of 1 electron and q quasiholes (or "vortices", or 
"flux quanta") at filling factor v — 1/q, (i) is a boson; (ii) sees no net effective magnetic field. 
When both (i) and (ii) are true, the composite can "Bose condense", i.e. have long-range 



order. For the Laughlin states, (i) and (ii) follow from the calculations in |15 



More formally, let if)* create an electron in the lowest Landau level (LLL), and let 



A' 



U(z) = U(z-z i ) 



(2.6) 



i=l 



be Laughlin's quasihole operator (acting in the LLL). If |0l; N) is the normalized Laughlin 



state for N particles, then one can show [12 



lim lim (0 L ; N\ U\z) q i)(z)^(z')U(z') q |0 L ; N) e -*W -*I*T = p. 

\z— z'\— +oo N^oo 



(2.7) 



Given our choice of order parameter operator, the Laughlin state is the state with the most 
order. In fact, we have [EJ 



|0 L ;AT) = (j d 2 z^{z)U(z) q e- 1 ^ 2 } |0> 



(2.8) 



where |0) is the vacuum (no electrons). This says that the Laughlin state is precisely a Bose 
condensation of the composite bosons. 

The order parameter in a general state obeys JCJ a system of Landau-Ginzburg-Chern- 
Simons equations [Hj], which also involve internal vector and scalar potentials. Similar 
equations are found in another, more field theoretic approach [n, |I71 which has been reviewed 
recently by Zhang The ideas have been extended [17, [TP], ^TJ to the hierarchy scheme 



and its generalizations [TP], which give states with abelian statistics for all 



filling factors p/q < 1. These also have interpretations as CFT correlators 0, TP ]. 



III. Pfaffian State and Nonabelions 
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Consider the wavefunction 



Vnizi, •■-,z N ) = Pfaff f— L_ ) TJ( Zi _ Zj )* e -\W (3.1) 



The Pfaffian is defined by 



I L/2 



PfaffMiJ = 2V*(L/2)\ ^ Sgna II M^i M 2k) (3.2) 

v / /' ctGSl k=l 

for an L x L antisymmetric matrix whose elements are My, or as the square root of the 
determinant of M; S n is the permutation group on n objects. \l/pf can be regarded as a 
wavefunction for spinless or spin polarized electrons in the LLL if q > is even, since then 
it is antisymmetric; the filling factor is v = 1/q. Note that Pfaff ((zj — Zj)~ r ), r odd, < q 
would also give a valid wavefunction. 

The idea behind the construction of this wavefunction was the following. From the 



calculations in fL5| , it follows that, in an incompressible fluid state of electrons at filling 
factor 1/q, the composite ^{z)U{z) q is a neutral boson if q is odd (as mentioned in §11) and 
is a neutral fermion if q is even. Hence, for q odd the operator can Bose condense, which 
gives the Laughlin states. For q even, on the other hand, it cannot condense singly, but 
pairs of fermions can condense, as in the BCS theory of superconductivity. For the spinless 
or spin-polarized case, the pairing function must be of odd parity to satisfy Fermi statistics. 
Thus a possible function is 

1 \ N / 2 

d 2 zd 2 w 4>Hz)U q (z)4>Hw)U^w)e--^ z \ 2+ \ w \ 2) ) 10). (3.3) 

z — w J 



Writing out this function in co-ordinate representation yields (3T). The Pfaffian is precisely 
the sum of products of fermion pairs, antisymmetrized over all distinct ways of pairing. More 
generally, the pairing function (z — w)' 1 can be replaced by {z — w)~ r with r odd, provided 
r < q, as noted above. This construction was also inspired by noticing that the Haldane- 



Rezayi (HR) spin-singlet state for electrons at v — 1/q, q even |25], can be written in the 



form involving det(zj — Zj)~ 2 - The real-space form of the spin singlet BCS wavefunction 
for spin-i fermions is just such a determinant of an even parity pairing function, as is well 
known (it is perhaps less well known that the analogous result for spinless or spin-polarized 
fermions is a Pfaffian). Thus the HR state is a condensate of pairs of spin-| composite 
neutral fermions, ipt( z )U(z) q ||! 

The pairing picture suggests two kinds of possible excitations. One is the analogue of the 
BCS quasiparticle, obtained by adding composite fermions, or by "breaking pairs". Thus 
a state with such a neutral fermion excitation localized at z is obtained by acting on the 
ground state with ijft (z)U(z) q . (In the HR state, one likewise obtains spin-| fermions.) 
These excitations have no analogue in the Laughlin states since there the bosons are already 



condensed; the only excitations, other than the collective density mode [26], are Laughlin's 
quasiparticles, which are fractionally charged. In the order parameter picture, these are vor- 
tices, i.e. the order parameter winds in phase by a multiple of 2ir around each quasiparticle, 
which thus resemble flux quanta in a superconductor, the flux quantum being $o = hc/e 
since the order parameter carries charge 1 from the single electron that it contains. The flux 
$ determines the charge e* through the quantized Hall relation, e* = i/($/$o). Similarly, in 
the paired states, one expects to find excitations corresponding to flux quantized in multi- 
ples of |$0; since the order parameter contains two electrons, and these will have charge in 
multiples of (2q)^ 1 . 

These ideas motivated the following wavefunction 0] for a pair of quasiholes in the 
Pfaffian state: 

^Pfaff+qholes^l, • • • , Zn] f 1, V 2 ) = (c 

J2 Sgna nS[(M2fc-i) - vi)(z a(2k) - v 2 ) + (vi g v 2 )\ I ^ )t?c „iy- ;kd 2_i M 2_i 
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The factor in brackets can be rewritten as 

( Vl - , 2 )^Pfaff f ^-^-^) + ^^^) 

(The reason for the assumed form of the exponent of {v% — v?) will be explained below.) The 
extra factor in each term of the Pfafnan resembles a pair of Laughlin quasihole operators, 
except that in each term each factor acts on only one member of each pair of fermions. 
Therefore, in an average sense, each quasihole is a half flux and carries charge l/2q. Note 
that this construction cannot produce one quasihole. One way to see that this is impossible 
is that working on a compact geometry like the sphere, the total flux is quantized in units 
of $o- Another feature of ( |3.5|) is that as v\ — > t>2, we recover a Laughlin quasihole of charge 
1/q. 

While our construction of the Pfafnan state was motivated by order parameter consider- 
ations, we can also give an interpretation using conformal field theory, which then suggests 
the full structure of the system of excited quasihole states. Introduce free, massless real 
(Majorana) fermions x in 1 + 1 dimensions: 

<X(*)XC0> = —, (3-6) 
z — z 

Then using also the free scalar field as before, construct 




This reproduces the FQHE wavefunction ( |3.1|) . To reproduce the two quasihole function 
( p.5|) , we should understand the CFT we are dealing with. Majorana fermions arise naturally 
in the two-dimensional Ising model at its critical point. The fermions, in conjunction with 
their leftmoving partners represent the energy density fluctuation e(z,z) = x(z)x( z ) 

which hence has dimension x — 1 and so the correlation length exponent is v — (2 — x)^ 1 = 1. 
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The other fields in the Ising model are the spin fields, i.e. the Ising spin itself. These are 
operators which produce a square root branch point in the fermi field 0. With the latter 
we can reproduce the two quasihole state: 

Wholes = (a( Vl )e^ M a(v 2 )e^ M f[ x (z i )eW«) e ^S*''*«'^ . (3.8) 

The square roots produced by the spin fields a are cancelled by those produced by the "half 
fluxes" exp(-^=(p(vi)), so the electron wavefunctions are single-valued. The factor {v\—v 2 )~^ 
in (|3.5|) is now explained; it was inserted to make the equality with the correlator hold, and 
reflects the fact that the conformal weight of the spin field is 1/16 (and hence rj = 1/4 in 
the Ising model). 

The generalization to 2n of our "half flux" quasiholes now seems self-evident: we should 
insert 2n of the combinations aexp(^=if) into the correlator. However, it is known 
that the Ising model part of such an expression is ambiguous; for fixed positions v i, ... v 2n 
of the spin fields there are many different possible correlators (conformal blocks), forming 
a vector space of dimension 2 n_1 (over the complex numbers) independent of how many 
fermions x are present. To resolve the ambiguity of notation, the spin fields should be 
replaced by "chiral vertex operators" J8]. Furthermore, the monodromy of these functions 
is nonabelian, i.e. analytically continuing Vi around Vj produces a linear combination of 
the original branches of the functions, the coefficients being some braiding matrices that 
don't commute. At present we do not know the explicit functions for 2n spin fields and N 
fermions in general, except in the case N = 0. The case N = may seem strange from the 
point of view of electron wavefunctions, but it does at least provide an explicit realization 
of nonabelian statistics, which is just that given in §1, where the wavefunctions for the "ex 
particles" are just N = Ising conformal blocks. In fact, the braiding properties of the spin 
fields are independent of N. Note that, when discussing the N fermion functions as electron 
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wavef unctions, we must show that they are 2™ _1 linearly independent functions of the electron 
co-ordinates but this is easily done using the operator product expansion of the spin fields 
as the co-ordinates v±, ... v 2n approach each other in pairs V\ — > v 2 , v 3 — > t> 4 , etc. [^TJ, even 
without full knowledge of these functions. The same operator product expansion also shows 
that as any Vi approaches any Vj, the leading term is equivalent to a Laughlin quasihole of 
charge 1/q, which is clearly a desirable property. This completes our construction, though it 
remains to check that adiabatic transport of our quasiholes does give the same nonabelian 
statistics as the monodromy of the conformal blocks. We are confident that this is true 
because of the existence of 2 n_1 functions and because analogous properties hold in the 
abelian examples of §11. 



Before closing, some comments on recent related work. Wen [28 has argued that exci- 



tations of wavefunctions arising in one of Jain's later constructions |24j possess nonabelion 
excitations, using the point of view of ref |§ of wavefunctions as conformal blocks. These 



examples involve SU(N) symmetry. Greiter, Wen and Wilczek |29| have studied numerically 
the vicinity of v = 1/2 for certain Hamiltonians in spherical geometry, and find evidence for 
an incompressible state with the excitations (neutral fermions and charge 1/4 quasiholes) of 
the Pfaflian state. However, they claim that the quasiholes obey abelian 6 /it = 1/8 statistics. 
Their arguments are based on an "adiabatic heuristic" principle that connects the FQHE 
state to a paired state of fermions in zero magnetic field given by the Pfaflian alone. While 
this picture is extremely similar to the order parameter or condensation picture put forward 



in previous papers p2| , |13| , [14] , [| and herein, we do not agree with the conclusion that the 
statistics is abelian. The full structure of the Landau-Ginzburg-Chern-Simons theory of the 
state must account for the fermion excitations, and not only the paired order parameter. 
Their conclusion is what one would obtain by neglecting the Pfaflian completely when cal- 
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culating the statistics of the quasiholes by adiabatic transport a la [To]. In fact, while a 
direct search for nonabelian statistics must use wavefunctions for 4 quasiholes (which are 
not given by Greiter et al), there is already a difference in the apparent abelian statistics for 
two quasiholes, as shown by our eq. ( |3.5|) , which implies 9/n = 1/8 — 1/8 = 0; this effect 
would be due to the Pfaffian. These authors also point out that the Pfaffian state is the 
unique highest density zero energy eigenstate of a certain three-body Hamiltonian. 

Clearly, once the possibility of nonabelions is accepted, many questions remain. How 
would we confirm their statistics in an experimental setting? Is a gas of nonabelions a 
superfluid, as for anyons fl30f ? 
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